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In an experiment in the Institute of Continuous Media Mechanics in Perm (Russia) an non-stationary screw dynamo is intended to be 
i; realized with a helical flow of liquid sodium in a torus. The flow is necessarily turbulent, that is, may be considered as a mean flow and 
a superimposed turbulence. In this paper the induction processes of the turbulence are investigated within the framework of mean— fleld 
I I electrodynamics. They imply of course a part which leads to an enhanced dissipation of the mean magnetic field. As a consequence 
^ of the helical mean flow there are also helical structures in the turbulence. They lead to some kind of o-effect, which might basically 
support the screw dynamo. The peculiarity of this Q-effect explains measurements made at a smaller version of the device envisaged 
^ for the dynamo experiment. The helical structures of the turbulence lead also to other effects, which in combination with a rotational 
I J shear are potentially capable of dynamo action. A part of them can basically support the screw dynamo. Under the conditions of the 
^ , experiment all induction effects of the turbulence prove to be rather weak in comparison to that of the main flow. Numerical solutions of 
i: the mean— fleld induction equation show that all the induction effects of the turbulence together let the screw dynamo threshold slightly, 
• ' at most by one per cent, rise. The numerical results give also some insights into the action of the individual induction effects of the 
turbulence. 
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A screw-like, that is, helical motion of an electrical l y con ducting fluid is capable of dynamo action. This 
has been first theoretically shown by Ponomarenkol (|l97,^ ). Later, in 1999, it was ex perimentallv demon- 



OO \ strated i n the liq uid sodium facility of the Institute for Physics in Riga (Latvia) (|Gailitis et al\ 20001 . 



O . 12001 al lhl. l2002al lbL 1200.1 l2004l V It is also the background of another liquid sodium experiment which is 
. under preparation in the Institute for Contin uo us Media Mechanics in Perm (Russia) (Denisov et al. 1999, 
\0 Frick et al. 2001, 2002) ()DeTiisov et »i.lll 9991 . IPHck et rfi.ll200l1 . l2002h . 



In this experiment a torus-shaped vessel filled with liquid sodium is rotated about its symmetry axis 
Yl , and then suddenly stopped. After the stop a flow of sodium occurs inside the vessel and decays in the 
course of time. Due to diverters, that is, proper arrangements of blades inside the vessel, the flow becomes 
helical. In this way the possibility of a dynamo of Ponomarenko type arises for a certain time interval. 
r^ '' The large radius of the torus is about 0.40 m and the small one about 0.12 m. Rotation rates up to 50 rps 
Qh' are envisaged. Estimates show that magnetic Reynolds numbers sufficient for dynamo action are indeed 
J> \ feasible. 

The critical magnetic Reynolds number is definitely higher than 10. For liquid sodium the hydrodynamic 
Reynolds number exceeds the magnetic Reynolds number by a factor of about 10^. Therefore in the case 
of a dynamo the fiuid flow is necessarily turbulent, that is, may be considered as consisting of a main flow 
and a superimposed turbulence. The turbulence affects the dynamo in two ways. Firstly it influences the 
profile of the main fiow. Secondly it leads to small-scale induction effects. In particular the turbulence 
will enlarge the effective magnetic diffusivity of the fiuid, what raises the dynamo threshold. Likewise 
other small-scale induction effects will occur. As a consequence of helical features of the turbulence an 
a-effect or related effects are possible. Basically they could also support the screw dynamo or even open 
the possibility of another type of dynamo. 
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In this paper we want to estimate these direct influences of turbulence on the electromagnetic field on the 
basis of mean-field elect rodynamics. We will rely on a recent general calculation of the mean electromotiye 
force due to turbulence ( Eadler and Stepanov 200.^ . which was partially motivated by this experiment. In 
section 121 we describe the general concept of our investigation. In section El we discuss possible structures 
of the mean electromotive force which can be concluded from symmetry arguments, and in section |31 we 
explain results of more detailed calculations. On this basis we deliver in sectionOestimates for the influence 
of the turbulence on the dynamo. Finally in section El the main results are summarized. 



2 The mean— field concept 

Let us first recall some ideas of the mean-field conc ept which are important for the following (see, e.g.. 



Krause and Riidleilll97ll . iMoffattllToT^ . iRiidleil 1200(1 ). We start from the induction equation that governs 



the behavior of the magnetic field B in an electrically conducting fluid, 

- V X (U X B) -r/V^B = 0, V-B = 0. (1) 

Here U is the velocity of the fluid and ij its magnetic diffusivity, which is assumed to be constant. 

The velocity U and, as a consequence, also the magnetic field B are assumed to consist of large-scale 
parts and small-scale turbulent fluctuations. We deflne the mean magnetic and velocity flelds B and U as 
averages over space or time scales larger than those of the turbulence and put B = B + b analogously to 
U = U + u. We assume that the Reynolds averaging rules apply. Taking the average of equation we 
obtain the mean-fleld induction equation 

9iB - V X (U X B + £) - T]V^B = , V • B = , (2) 

where £ is an electromotive force due to the fluctuations of magnetic fleld and velocity. 



f: = uxb. (3) 

The equation for b resulting from and ((SJ allows us to conclude that £ is determined by U, u and B. 
More precisely, £ at a given point in space and time depends not only on U, u and B at this point but also 
on their behavior in a certain neighborhood of this point. We adopt the frequently used assumption that 
B varies only weakly in space and time so that £ in a given point depends on B only via its components 
and their first spatial derivatives in this point. Then £ can be represented in the form 

Si = aij Bj + bijk dBj/dxk , (4) 

with tensors aij and hijk being averaged quantities determined by U and u. Here and in the following we 
refer to a Cartesian co-ordinate system {xi,X2,x^) and adopt the summation convention. Relation (jlj is 
equivalent to 

£ = —a o B — 7 X B 



— Q o B — 7 X B 
-/3 o (V X B) - d X (V X B) - K o (VB)(^) 



(5) 



(see Radler 1980, 2000), where a and (3 are symmetric tensors of the second rank, 7 and 5 are vectors, 
and K is a tensor of the third rank, all depending on U and u only. (VB)'^'*) is the symmetric part of the 

gradient tensor of B, i.e. (VB)j^^ = ^{dBj/dxk + dBk/dxj). Notations like a o B or k o (VB)^'^^ are used 

in the sense of {a o B) j = ctijBj or (k o (VB)(^))i = Kijfe((VB)(^))jfc. The term with a in Q describes the 
a-effect, which is in general anisotropic, that with 7 a transport of mean magnetic flux by the turbulence. 
The terms with /3 and S can be interpreted by introducing a modifled magnetic diffusivity, again in general 
anisotropic. They are usually accompanied by the k term, which allows no simple interpretation. 
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3 The structure of the mean electromotive force £ in the presence of a screw motion 
3.1 Mean motion 

We consider the Ponomarenko dynamo here in its original cylindrical geometry. As observed in other 
investigations this provides a reaso nable appr oximation for the situation in a torus even for aspect ratios 
as large as in the planned device (|Stena,novir2nnnl . We assume that it is primarily a helical mean flow 
described by the mean velocity U which generates or maintains the mean magnetic field B. Referring to 
a proper cylindrical co-ordinate system (r, (/?, z) we put 

U= (0,I7^(r),I7,(r)), U^ = ^{r)r. (6) 

As mentioned above we will study the influence of the turbulence, that is, of the velocity field u, on the 
mean magnetic field B in the framework of mean-field electrodynamics. 

For a more detailed discussion of the mean electromotive force £ caused by the turbulence and for its 
calculation we change temporarily from the frame of reference in which ® applies to another one rotating 
about the axis r = and moving along it. Let us consider £ in a given point at the surface r = rg. 
The moving frame is then fixed such that the mean velocity of the fluid at the surface r = rg is equal to 
zero. The result for S found in this rotating system, considered as a vector with the usual transformation 
properties, applies in the original frame of reference as well. In the moving frame the mean velocity U has, 
again in the cylindrical co-ordinate system introduced above, the form 

U = (0, (17(r) - f](ro))r,I7,(r) - U M)) ■ (7) 

In contrast to the original frame, in the rotating one a Coriolis force occurs, which is defined by an angular 
velocity fi, given by 

n = n{ro) z , (8) 
where z means the unit vector in the z direction of the cylindrical co-ordinate system. 



3.2 Homogeneous turbulence 

Consider now € in the moving frame at a point with U = 0. Then £ and so the quantities Uij and bij^, 
or OL, 7, /3, S and k, depend on U and u in a certain neighborhood of this point. The dependence on u 
implies again a dependence on U and of course also a dependence on ft. We restrict our attention now to 
the case in which U varies only weakly in the relevant neighborhood and assume that its behavior can be 
there sufficiently precisely described by the gradient tensor VU, with respect to Cartesian co-ordinates 
defined by (VU)ij = dUi/dxj. In addition we assume until further notice that the turbulent fluctuations u 
deviate from a homogeneous isotropic mirror-symmetric turbulence only as a consequence of the gradient 
VU of U and the Coriolis force defined by fl. We may split VU into its symmetric and antisymmetric 
parts. The first one is the rate of strain tensor D, that is, Dij = ^{dUi/dxj + dUj/dxi). It describes the 
deforming motion near the point considered. Due to the incompressibility of the fluid we have V • U = 
and therefore Da = 0. The second part, A, given by Aij = ^{dUi/dxj — dUj/dxi), corresponds to a rigid 
body rotation of the fluid near this point. We may represent it according to Aij = —^eijiWi by a vector 
W = V X U. Until further notice we do not use specifications of D and W according to 0. They will be 
introduced only later. 

We now utilize the symmetry properties of the basic equations governing B and U, or B, U, b and 
u, from which S and so the quantities a p and buk, or ct, 7, , S and k, have to be derived. They allow 
two important conclusions (see also, e.g., Eadler and Stepanovli2nfl 5). The first reads that these quantities 



cannot contain any other construction elements than the isotropic tensors 6im and eimn, the vectors fl 
and W and the tensor D. The second uses the distinction between polar and axial vectors as U and B, 
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respectively, and between true and pseudo tensors. It reads that Uij and bijk, further a and k are pseudo 
tensors, /9 is a true tensor, 7 a polar vector, and S an axial vector. It is then important that fl and W 
are axial vectors and D is a true tensor. 

Let us consider first aij and 7^. The mentioned construction elements 6im, ^imn-, ^1 W and D allow us 
neither to built a pseudo tensor of the second rank nor a polar vector. That is, we have 



a,; 



, 7* = . 



(9) 



We stress that this applies independently of any assumptions concerning linearity in ^l, W and D like 
those we will introduce below. 

In contrast to this there are several non-zero contributions to (3ij, 5i and Kijk- For the sake of simplicity 
we give here only those of them which are linear in CI, W and D, 

i^ijk = iK^^'^njd.k + ^Aj) + \n^^\Wj5ik + Wk6,j) + K^^'^e.jiDki + eikiD.i) . (10) 

The coefficients 6^^\ ■ ■ ■ K^^^ are determined by u but do not depend on fi, W and D. Because 

of V • B = 0, terms of Kijk containing 6jk would not contribute to £ and have therefore been dropped. 
As a consequence of Q and (|1()|) we have 



-/3(°) V X B - /3(^)d o (V X B) 
-{6^^^n + 6^^^W) X (V X B) 

-(K(^)n + K^^^W) o (VB)(^) - k(^) k(D) o (VB)W , 



:ii) 



where k(D) is a tensor of the third rank defined by kijk = eijiDik + eikiDij. 

The and /?(^) terms in (jllj) make that the mean-field diffusivity deviates from the original, that 
is, molecular magnetic diffusivity of the fluid. Due to the term the mean-field diffusivity becomes 
anisotropic. The (5*^^^ and (5^^) terms, too, can be discussed as contributions to the mean-field diffusivity. 
They lead to skew-symmetric contributions to the diffusivity tensor. In another context the effect described 
by the (J*-^) term has been called "fi x J-effect". It has been shown that this effec t in combination with 
a differen tial rotation, her e a depende nce of Q or r , is able to establish a dynamo ( Radlei 19691 . Robert d 
IT972 . .Moffatt and Procton 1982. Rii dler et al\\2m±. The 5^^^ term describes an effect analogous to the 
n X J -effect, which has been revealed only recently ( Roeachevskii and Kleeorinll2003l . Radler and Steoanov 
200i). It occurs even in the absence of the Coriolis force as consequence of a shear in the mean motion. The 



possibility of a dynamo due to this effect is still under debate, see section 15.31 We refrain from discussing 
the K^^\ K^^^ and n^^^ terms here. 



3.3 Inhomogeneous turbulence 

Let us now relax the assumption that the original turbulence is homogeneous and isotropic. We admit an 
inhomogeneity and an anisotropy due to a gradient of the turbulence intensity and introduce a vector g 
in the direction of this gradient by putting Vu-^ = gu^. Then we have to add g, which is a polar vector, 
to the above-mentioned construction elements of Oij, 7^, Pij, 6i and Kijk- As a consequence cx and 7 can 
well be non-zero. For the sake of simplicity we assume that the influence of g on a, /3, 7, S and k is so 
weak that they are at most of first order in g. We have then 

aij = af\g ■ n)6^j + af\ginj + g^Qi) + a^^\g ■ W)5,, + a^^\giWj + g^Wi) 

{^iimDji + ejimDii) gm (12) 
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whereas (jlUp remains unchanged. 
The general form of £ reads now 

S = -a?^ (g • n)B - af^ ((n • B)g + (g • B)^) 

-af'\^ ■ W)B - 4^^((W • B)g + (g • B)W) 
-a(^)Q(g,D) oB 

_(^{o)g + ^(n)g ^(H/)g X W + 7(^)g o D) X B 

V X B - /3(^)d o (V X B) 
-(<5(^)n + ,5(^)w) X (V X B) (13) 
-{K^^^fl + k('^)w) o (VB)(^) - k(^) k(D) o (VB)W , 

where a{g, D) is a symmetric tensor defined by Oij = (eumDij +ejimDii)gm- In the fohowing the induction 
effects described by the individual contributions to € are called a*^^^-effect, a^^^-effect, • • • K^^^-effect or, 
more summarizing, a-effects, /^-effects, • ■ • /t-effects. 

The occurrence of an a-effect as a consequence of the gradient of the turbulence intensity and the Coriolis 
force, that is, of g and ft, is known for a long time. The possibility of an g-effect due to a combination 
of g and W, or of g and D , has however been revealed only recently ()Rogachevskii and KleeorinlEnn^ . 



Radler and Stepanov 200^. Likewise the 7-effect and its modification by the Coriolis force, that is fl 



is known for a long time, but modifications by W and D have been found only in the last-mentioned 
investigations. 



3.4 a effect 

We still remain in the co-moving frame defined above, in which U is given by and consider a point 
with r = tq. As a consequence of ((TJ both W and D take specific forms. The non-zero components of 
W are = —dUz/dr and Wz = rdO/dr. As can be concluded with the help of the relations in the 
Appendix the only non-zero elements of D are Dr^p = D^r = § dQ/dr and Drz = Dzr = \ dU z/dr. As for 
n we refer to (jHJ. 

In view of the Perm dynamo experiment the simplest assumption on g is that it is a radial vector. In this 
case the tensor cx has a remarkable peculiarity. The terms with af^^ and a[^^ disappear. In the cylindrical 
co-ordinate system cx has the form 



(0 oi 
ai 
02 / 



ai = -(«r-^a(^))^^^, a2 = (ar-^a(^)).^ + af (14) 

where U^p and U z are, of course, measured in the rotating frame of reference. All diagonal elements of o; 
are equal to zero. Therefore a mean magnetic field in z direction does not generate any mean electromotive 
force in (/? or z direction. The last statement remains true even if we relax the assumption on linearity with 
respect to g. 

In this context experiments with torus-shaped devices similar to but smaller than that prepared for 
the dynamo experiment are of interest. In an experiment with water instead of sodium clearly helical 
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small-scale motions have been observed (jFrick et »i.ll2nn2h . This suggests the existence of an a-effect in 
the conducting fluid. In the case of a complete a-effect an azimuthal mean magnetic field leads to a mean 
electromotive force with a non-zero azimuthal component. However, our above result concerning a, with 
the z direction interpreted as the azimuthal one, does not predict such a component. The existence of 
helical motions alone is not sufficient for the occurre nce of a cornplet e a-effect. Another experiment with 
gallium instead of sodium was carried out, too (Nos kov et all It clearly confirmed that there is 

indeed no azimuthal electromotive force or electric current due to an azimuthal magnetic field. 



4 Specific results for the mean electromotive force £ 

Our considerations on the structure of the electromotive force £ did not provide us with relations between 
the coefficients af^\ a'^\ k^'^'^ and the parameters of the turbulent flow. We r e lv he re on recent 



calculations of these coefficients presented in detail in a paper by iRadler and Stepanovl ()2nnfil ) . They were 
carried out just under the restrictions introduced above, namely that £ is only weakly influenced by the 
Coriolis force and the gradient of the mean velocity so that it is linear in ft, W and D. In addition, again 
in the sense of restrictions already made, it was assumed that all mean quantities determining £ vary only 
weakly in space and time so that it is also linear in the operator V acting on quantities like or B, and 
contains no time derivatives of them. Finally the second-order correlation approximation was used. In this 
approximation the velocity correlation tensor Vij of the second rank. 



Vijin, T; r, t) =Ui(R + r/2, T + t/2)uj(R - r/2, T - t/2) , (15) 

plays a central part. The calculations have partially been carried out in the Fourier space, and the Fourier- 
transformed correlation tensor Vij has been used, deflned such that 

tijj (R, T; r, t) = jj i;jj(R, T; k, exp(i(k • r) — Lijt) d'^/c dw . (16) 

The "original" turbulence, which occurs in the limit of vanishing Coriolis force and vanishing gradient of 
the mean motion, is assumed to deviate from a homogeneous isotropic and mirror-symmetric one only by 
a gradient of its intensity. The corresponding Fourier-transformed correlation tensor Vij was assumed to 
have the form 

i}i,-(R,r;k,c^) = ]^{5ij - ^ + l^ih^j - kjVi))W{Yi,T-k,uo) , (17) 

with some function W describing the kinetic energy distribution in the Fourier space. It was chosen such 
that 



u(R + r/2, T + t/2) ■ u(R - r/2, T - t/2) = u2(R, T) exp(-r72A^ - t/|Tc|) , (18) 

where v? describes the turbulence intensity in the limit of vanishing Coriolis force and mean velocity 
gradient, and Ac and Tc are correlation length and time in this limit. 

The results for the coefficients a^^\ ck2^\ • ■ k^^-* can be represented with the help of the Prandtl 
number Pm of the fluid, the magnetic Reynolds number i?„ of the turbulent motion and a dimensionless 
parameter q, 

Pm = v/'q, Ru = uXc/r], q = Xl/r]Tc, (19) 

where u = \/u^. Clearly Ru, like u, in general depends on position. The quantity q is the ratio of the 
magnetic diffusion time Xc/v to the correlation time Tc- We speak of low-conductivity limit ii q —^ 0, and 
of high-conductivity limit if g ^ oo. 
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Table 1. Numerical values of Q^^\ ^2^'i ' ' 'j ^^^^ for various values of P^ti and q. 



p 


1 


10-i 




10^'' 


10~* 


10~^ 


10~^ 


10^^ 


10^'' 


q 


0.02 


0.02 


0.02 


0.02 


0.02 


0.002 


0.02 


0.2 


0.02 




n n7Q9 


u.Doy 




i .0 




u.yoz 


y.ui 


o I . 1 


y.uo 


-In) 


-0.104 


-0.654 


-1.98 


-2.35 


-2.2 


-0.195 


-2.18 


-22.5 


-2.17 


(W) 

a] ' 


0.0439 


0.38 


1.82 


3.76 


4.38 


0.46 


4.46 


43.6 


4.47 




-0.0094 


-0.139 


-0.57 


-0.972 


-1.07 


-0.092 


-1.08 


-10.9 


-1.08 




-0.0437 


-0.171 


-0.374 


-0.191 


-0.048 


-0.015 


-0.027 


-0.266 


-0.025 




0.0546 


0.0546 


0.0546 


0.0546 


0.0546 


0.049 


0.0546 


0.0555 


0.0546 




-0.004 


-0.028 


-0.067 


-0.089 


-0.093 


-0.066 


-0.093 


-0.105 


-0.094 




-0.0086 


-0.0541 


-0.409 


-1.12 


-1.39 


-0.165 


-1.42 


-13.5 


-1.43 




-0.0888 


-0.766 


-3.09 


-5.44 


-6.05 


-0.629 


-6.13 


-60.3 


-6.13 


,3(0) 


0.109 


0.109 


0.109 


0.109 


0.109 


0.0979 


0.109 


0.111 


0.109 


j3(.D) 


0.059 


0.558 


2.03 


3.09 


3.25 


0.293 


3.27 


33.1 


3.27 




-0.004 


-0.028 


-0.067 


-0.089 


-0.093 


-0.066 


-0.093 


-0.105 


-0.094 




0.0214 


0.0214 


0.0214 


0.0214 


0.0214 


0.0125 


0.0214 


0.0255 


0.0214 




0.0072 


0.098 


0.268 


0.364 


0.382 


0.24 


0.386 


0.25 


0.386 




-0.0052 


0.0468 


0.14 


0.19 


0.202 


0.13 


0.204 


0.236 


0.204 




0.114 


0.631 


2.13 


3.21 


3.38 


0.373 


3.39 


33.2 


3.4 



As for the general results concerning af^ , 02 ' ' ' '^^^^ refer to the mentioned paper, in particular 
relations (55)-(58) and figures 1 and 2. These results cover arbitrary Pm and q. We may write af^\ a'^\ 
• ■ •, K^^^ except 7'''^'* and Z?'-'^-* in the form 

af ) = a?\Pr^, q) RI \l , = af (P^, g) i?^ , • • • , k^^) = k^^^P^, q) RI Xl , (20) 

and 

1^'^=1^'\q)RIv, P^^=i3^'\q)Rlv. (21) 

The a^^\ (x^\ ■ • •, kS^^ are dimensionless functions of Pm and q. Only 7^*^^ ^^^^ and 5^^) are independent 
of Pm- In the case of 7^^^^ and (3'^'^'^ this independence appears as natural when considering the derivation 
of these results, in the case of J(^) it occurs as a accidental compensation of different influences. Further 
we have 7^'''' = ^ j3^^'^ and t'*'^^ = (J*-^-* . As will become clear in the following section here only results for 

specific Pm and q are of interest. Table 1 gives some values of the a^^\ ct2'\ • ■ ■, R^^^ for such Pm and q. 



5 The effect of turbulence in the Perm dynamo device 
5.1 Estimates of the turbulence effects 

Let us make a few estimates of the turbulence effects to be expected in the Perm dynamo device on 
the basis of the general form of the mean electromotive force f , the relations (|2U|) and (|21|) . and 
numerical values represented in table 1. As for the properties of the liquid sodium we put rj = O.lm^/s 
and Pm = 10~^. As a typical value U of the mean velocity U we use U = 50m/s. For the velocity u of the 
turbulent motion, the correlation length Ac and the correlation time Tc we assume u = Im/s, Ac = 0.01m 
and Tc = 0.05 s. This implies Ru = 0.1 and q = 0.02. Further |g| is estimated by 1/i?, where R means the 
radius of the fiuid cylinder. Likewise |ri|, |W| and |D| are estimated by U/R, and spatial derivatives of B 
by B/i?. We put R = 0.12m. _ 

Relation ((Ti^ for £ ignores any time variation of B. Since we are in the limit of small q this is justified 
as long as the characteristic time of this variation is large compared to A^/ry, with the above data 10~'^s. 
The calculations described below show that this characteristic time is not smaller than 10~^ so that this 
condition is satisfied. 

Consider first the /3-effects, which may be interpreted by introducing a mean-field version of the magnetic 
diffusivity instead of the molecular diffusivity i]. This mean-field version is in general a tensor, r^mij = 
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(r/ + /3(°))(5ij + The contribution of the /3(°)-effect to the mean-field diffusivity, measured in 

units of T/, is f3^^^ /r] = j3^^^R^. With the above data it takes a value of about 10~^. The corresponding 
contribution of the /S^^^-effect is of the order of /S^^^U/rjR = (3^^\UR/r]) {X^/RfRl, with the above data 
about 1.3 • 10-2. 

It is natural to compare the 7-effects, which describe a transport of mean magnetic flux, with that of the 
mean velocity U. As for the 7''°^-effect we find that |7'''''*g|/lU| is of the order oij^^^/UR = ^^^\ri/UR) i?^, 
which takes a value of about 10"^. In analogous estimates for the 7^^^, 7^^^ and 7^^-' -effects the ratios 
corresponding to j^^^/UR are 7^^^/i?^, 'j^^^/R'^ and 7^^^/i?^, all being independent of U. The first two 
are smaller than the last one, which is about 4 ■ 10"^. 

We may also compare the a, 5 and K-effects with that of the mean velocity U. They are then characterized 
by the ratios af^^ /R"^, a^'''^ /R^, a{^'' /R^, ■ ■ ■, 6'^^^/R'^, ■ ■ •, k^^^ /R'^, ■ ■ •, again being independent of U. 
As for the ai'^ /R?, /R^, ■ ■ ■ a^^^ /R?, the first one is about 6 • 10~^ but all others are smaller. We 
recall that the a\^^ or the a^'^^-effects vanish if g is orthogonal to or W, respectively. Further 6^^^ /R'^ 
is about 7 • 10"*^ but 6^^^ /R^ smaller. Among the k^^\ ■ ■ ■ k^-^) /R^ the last one is about 2 • 10""^ but the 
others are smaller. 

We may conclude from these estimates that the influence of turbulence on the screw dynamo in the 
Perm experiment should be rather small. There is hardly the possibility of another type of dynamo, which 
would become possible by turbulence effects. 



5.2 Dynamo model 

Some more detailed numerical investigations concerning the effects of turbulence on the screw dynamo 
have been carried out, mainly with the intention to find out in which sense they influence its excitation 
condition. Again g was considered as a radial vector such that Vu-^ = gn^. 

Let us first give equation governing B inside the cylinder and the supplementing relations in a 
non-dimensional form. In that sense we put 

n = non, u, = UoU,, u^ = ulx. (22) 

Here and Uq are constants describing typical values of O and Uz, e.g., being their values at r = 0, and 

and U z are dimensionless functions of r. Further uq means the value of u at r = 0, and x is again a 
dimensionless function of r. In the following we measure all lengths in units of the radius R of the cylinder 
and the time in units of R? jr]. Then ((21) can be written in the form 

9fB-Vx (i?c/UxB + i?2o(F(B)+^iit/G(B)))- v2b = 0, V-B = (23) 

with the parameters 

Ru = ^J^lR^ + UiRh, Ruo = uoXJrj, i = {\jRf. (24) 

and 

\J={QXn^{r)rXuUz) (25) 

F = -7(0)^1- X B-/3(°)xV X B (26) 

G = -af^ x'Cn^(Bzr + Brz) 

-(af ) - a(^)) x' iCn rQ' (B,r + 5,£) - (u U'z {B^i + Br(p)) 
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Figure 1. Profiles of Q, and (left), and of x (middle) and x' (right) for p = 1 (solid), p = 18 (dashed) and p = 100 (dotted) 

-P'-''^ iCn r n' ((V X B)^f + (V x B),(^) + (u K ((V x B),r + (V x B),z)) 

xCn^^ + S^"^^ xiCnrn'z- (u U', 0)) x (V x B) (27) 
xCn^i^- k^"^^ xiCnrn'i- Cu Ui<^))o ( VB)(^) 
X {Cn r n' ((VB) Wf - ( VB)(f)<^ _ (( VB)(^; - ( VB) g )z) 

-Cu U'. ((VB)Wf - ((VB)W - (VB)W)^, - (VB)gz)) 



Cn = 1/ + e , Cu = C/Vl + e , C = Uo/noR. (28) 

The r, (p and z are unit vectors in r, ip and z direction, and primes indicate differentiations with respect 
to the dimensionless radial co-ordinate r. Note that the 7'-''^, f3^^\ ••• k^^^ depend on q, most of 

them also on Pni. As for the components of (VB)^'^) in the cylindrical co-ordinate system we refer to the 
Appendix. Of course, equation has to be completed by boundary conditions. 
We will look for solutions of the above equations for B having the form 

B = B exp(i(m(/9 + kz) + Xt) (29) 



with a positive integer m, a real non-zero wavenumber k and a complex growth rate A. Note that by 
Cowling's theorem solutions with m = have to be excluded as long as Ruo = 0. 
As for Q and Uz we put 

0(.) = ^.(r) = ^°^^^^|^-^°^y-^ (30) 
cosh(/7j — 1 

with some parameter p. Interpreted in view of Uz, for p = 1 the flow has a Poiseuille profile, for p — > cx) 
a piston profile. The outcomes of the water experiments made to simulate the flow in the dynamo device 
fit best to p = 18. In view of y, which i s defin ed by ()22|) . we adopt a relation for the r.m.s. value u of the 
turbulent velocity given by Schlichtin^ (^6^), 



uE 

u{2.5ln— + 1.75) = Uz. (31) 

On this basis x has been calculated with u such that Pm = 10^^. Profiles of Cl and Uz and of x and x' 
for several p are depicted in figure ^ 

All calculations reported in the following have been carried out with Pm = 10^"' and q = 0.02. The 
corresponding numerical values of a2^\ a2^\ ■ ■ •, k^^^ are given in table 1. 

As for the surroundings of the fiuid cylinder two cases are considered. In the first one all surroundings 
are assumed to be free space. Since only non-zero k are admitted the total electric current through a 
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Figure 2. and k* versus f for the cases d = 0, p = 1 (dashed), d = 0.15, s = 5, p = 1 (dotted), and d = 0.15, s = 5, p = 18 (solid) 

in the absence of turbulence, that is, RuO = 

cross-section of the cylinder is equal to zero. Therefore it is to be required that B inside the cylinder 
continues as a single-valued potential field in outer space. In the second case it is assumed that the fluid 
cylinder is muffled by a rigid electrically conducting shell, in the experimental device made of copper, and 
there is free space outside this shell. Again B has to continue as a single-valued potential field in this free 
space. The thickness of the shell, measured in units of R, is denoted by d, and the ratio of the electrical 
conductivity of the shell to that of the fluid by s. The first case may be considered as the special case 
d = of the second one. The experimental device corresponds to d = 0.15 and s = 5. 

The equations (|23() ~()27 |) together with the ansatz ()29() have been reduced to a system of ordinary 
differential equation which, together with the boundary conditions, defines an eigenvalue problem with 
the growth rate A as eigenvalue parameter. After discretization of the equations this problem has been 
solved with the QR algorithm. In the following we restrict our attention to solutions with 5R(A) = 0. The 
corresponding values of Ru are denoted as critical values. 




5.3 Influences of the turbulence on the excitation threshold 

Let us first give critical values of Rjj for the ideal situation without turbulence, that is, Ruo = 0. For any 
given ( and p they still depend on m and k. We restrict ourselves to m = 1. All solutions of the equations 
considered with the same parameters but other m have higher critical values of Rjj- We denote the minimal 
critical Rjj for given ^ and p but varying k by R^ and the corresponding A; by A;*. Figure |21 shows i?^ and 
k* as functions of C for several cases concerning d and p. It shows in particul ar that a highly conducting 
shell around the fluid cylinder m arkedly favors the screw dynamo; (see also lAvalos-Zuniga et al. 20031 . 
Avalos-Zuniga and PlunianI 2005^ . 



We admit now turbulence. Since in the experiment C, will probably be close to unity we put now through- 
out = 1. As figure 131 demonstrates the turbulence generally enlarges the dynamo threshold. Consider 
the most realistic case d = 0.15, s = 0.5 and p = 18. Assuming that Ruo is of the order of 10^^ and ^ of 
the order 10~^ we conclude that grows by at most 0.5 percent. We point out, however, that all our 
calculations have been done with = 10~^ and q = 0.02. The data of table 1 suggest that the influence 
of the turbulence will be markedly larger if larger values of q have to be assumed. 

Several investigations have been made in order to find out which specific induction effects of the turbu- 
lence such as a-effects, 7-effects etc. hamper the screw dynamo or possibly support it. In that sense only 
parts of the mean electromotive force S have been taken into account. Both cases with p = ^- and more 
realistic ones with p = 18 have been studied. In the latter one parts of the profiles of fi, Uz and x ^'^^ 
markedly steeper. 

Let us consider first the a-effects. Note that the a^^'* and a[^^-effects are absent since g is orthogonal 
to both ri and W. The influence of the 0^2^^ and a^^-* -effects on the screw dynamo depends crucially 
on p. In cases with p = 1 it was found that they act against the screw dynamo. Figure Et- shows results 
for a more realistic case with p = 18, in which these effects clearly support the screw dynamo. 
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Figure 4. Ky versus Rno for cases in which no other turbulent induction effects than the a-effects (left) or the 7— effects (right) are 
taken into account, with 1^ = 1, d = 0.15, s = 5, p = 18 and different J (indicated by different line styles as in figureEJ 



Consider next the 7~effects. In all investigated cases, both with p = \ and also with higher they 
hamper the screw dynamo. Figure [IJa illustrates this for the above case with p = 18. 

Proceeding now to the /3-effects we note first that the P^^^-eSect may be interpreted in the sense of 
the mean-field magnetic diffusivity rj + P^^^ different from the molecular diffusivity rj. Since P^^'^ is always 
positive it raises the threshold of the screw dynamo. The P^^^~eSect makes the mean-field magnetic 
diffusivity anisotropic. It is then described by the tensor (ry + 6ij + P^^^Dij. In all investigated cases, 
with p = 1 and also with higher p, the P^^^-eSect for not too small ^ always dominates the P^^^-eSect 
and clearly supports the screw dynamo. This implies of course that the mentioned tensor is not positive 
definite; otherwise it had to act against the screw dynamo. Results for the above case with p = 18 are 
shown in figure EK- 

Both the 5^^^ and and also the K^^^-eSect always support the screw dynamo. The k^^^ and 

K^^^-effects and also all K-effects together act against it. This is partially reflected in figures ISh andlHK 
applying again to the above case with p = 18. 

We recall here that (5^^^ and K^^^-effects together with a sufficiently strong rotational shear open the 
possibility of a mean-field dynamo different from the screw dynamo, often called ft x J dynamo (Radler 
1969, 1980, 1986, Roberts 1972, Moffatt and Proctor 1982). The positive influence of the 6^^^ and k'^^^- 
effects on the screw dynamo, demonstrated by flgureEh, can possibly be interpreted in the sense of this 
dynamo mechanism. 

In contrast to the combination of the 5^^^ and K^^^-effects that of the 5^^^ and K^^^-effects, as figureEK 
shows, does not support the scre w dynamo. For the W x J dynamo, which has recently been proposed 
( Rogachevskii and KleeorinlEnn^ . the k^^\ P^^^ and K^^^-effects are important. However, it turned 
out that at least the simple model in plane geometry, which was used to explain this proposal, fails to 
work as a dynamo as long as 6^^\ k^^\ P^^^ and k^^^ are determined in the second-order correlation 
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approximation ( Radler and Steoanov 20051 . iRiidiger and Kitchatinov 20061 ^ . In this context it is of interest 
that the combination of the 6^^\ f3^^^ and ac^^^ -effects, as figure EJa shows, does not support the 

screw dynamo. This is also in so far quite remarkable as we already know that the /S^^^-eSect supports it. 
It is the K,^^^ and K^^^-effects which have a strong opposite influence. 

Our results show that the total influence of the induction effects of the turbulence on the threshold of 
the screw dynamo, which we have estimated above, is weaker than the influences of some of the individual 
effects. Obviously a part of them compensates each other. 
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6 Summary 

The electromagnetic processes in Perm dynamo experiment are influenced by turbulent motions accompa- 
nying the main flow. Their influence on the dynamo process has been investigated within the framework 
of mean-field electrodynamics. The critical magnetic Reynolds number of the main flow, which defines 
the onset of the screw dynamo, grows with the magnetic Reynolds number R^q of the turbulence and 
the also with the parameter that is, the ratio of the decay time of a magnetic field extended over one 
correlation length to the correlation time. Under the assumptions on these quantities made above, which 
we consider as realistic, R^j changes hardly by more than 0.5 percent. If they arc by any reason larger, Ry 
becomes larger, too. 

The electromagnetic effects of the turbulence are analyzed in terms of a, 7, 6 and K-effects. Despite 
helical structures in the turbulence there is no complete a-effect. The mean electromotive force due to this 
a-effect has not always a component in the direction of the mean magnetic field. This fact explains the 
measurements made at an experimental setup simulating the dynamo device in a smaller scale. Despite 
this peculiarity the a-effect supports the screw dynamo slightly. The 7-eflFect acts always against the 
screw dynamo. The two contributions to the /3-eflFect, the /^^^^ and the /?(^)-effect, act in opposite ways. 
The /3(°)-eflFect leads to an mean-field diffusivity larger than the molecular one and so raises the screw 
dynamo threshold. The P^^^ effect makes the mean-field diffusivity anisotropic and strongly supports the 
screw dynamo. Likewise the 5-effects support the screw dynamo whereas the K-effects act against it. It 
is known that the 5^^^ and K^^^-effects in combination with a proper, sufficiently strong rotational shear 
may work as a dynamo, sometimes called O x J dynamo. They also support the screw dynamo. Analogous 
to the X J dynamo the possibility of a W x J dynamo has been discussed in the literature, for which 
the 5'^'^\ k'^'^\ /3^^^ and K^^^-effects are important, but the question of its existence remained open so 
far. Interestingly enough it was found that the combination of these effects does not support the screw 
dynamo. 
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Appendix 

The symmetric part ( VV)^*) of the gradient tensor of a vector field V, in a Cartesian co-ordinate system 
defined by (VV)g^ = \{dVi/dxj + dVj/dxi), has in a cyhndrical co-ordinate system the components 



(VV)W 


= ^ ( 

dr ' 


(VV)W 


= (vv)(t.) 


(VV)g 


= ( vv)W 



^ y )^ dip dr r ' 

2^dz^ dr'' ^^^>^^~ r^ d^^""' ^ ' 



